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Since the first molecular systematic studies on marine sponges in the 1980’s, many papers 
have been published about levels of allozyme divergence between conspecific and 
congeneric sponge populations. Those genetic studies have indicated that sponges are more 
divergent than other marine invertebrates, a fact that was attributed to the high levels of 
genetic variation and morphological conservativeness found in Porifera. However, an 
analysis of 55 interspecific and 87 intraspecific pairwise genetic identity (/) values indicates 
a more complex picture. This study found that the average of / over all interspecific 
comparisons (/=0.42) was not much smaller than that found among other marine 
invertebrates (J=0.54), and the frequency distribution of J, for intraspecific comparisons, 
appears to be bimodal. Some genera were consistently highly divergent (/<0.30; 
Cinachyrella, Oscarella, Cliona, Spirastrclla and Tethya), whereas others were within the 
normal range of gene divergence (0.40 « 7 < 0.80; Chondrosia, Suberites, Petrosia, Plakina 
and Phyllospongia). Furthermore, in the genera Axinella. Chondrilla and Clathrina, both 
low and high levels of intrageneric genetic differentiation were found (0.13 < 7< 0.82). This 
pattern may reflect a large variance in the evolutionary age of genera in sponges, with very 
large levels of intrageneric gene divergence for some. We conclude with two non-mutually 
exclusive scenarios: a) genetic identity levels are too variable among sponge species to be of 
any use to evaluate taxonomic rank above species, or b) the range of evolutionary divergence 
in some genera of sponges is so broad that they may need revision. CJ Porifcra, gene 
divergence, allozymes, heterozygosity, molecular systematics, larval dispersal. 
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For marine organisms genetic markers have 
bcen extremely useful both for estimating levels 
of gene flow in structured populations (Burton, 
1996), and for the detection of sibling species 
(Knowlton, 1993; Thorpe & Solé-Cava, 1994). 
Allozyme electrophoresis has become the 
method of choice for alpha (i.e. at thc species 
level) molecular systematics of marine 
organisms (Thorpe & Solé-Cava, 1994; 
Knowlton & Weigt, 1997). The main advantage 
of allozyme electrophoresis for taxonomic 
studies is that it represents an independent sct of 
characters for the detection of sibling species 
(Solé-Cava & Thorpe, 1987). Genetic markers 
such as allozymes are particularly powerful for 
alpha-taxonomy (Hillis etal., 1996) because they 
can be used to detect reproductive isolation in 
sympatry (i.e. the biological species concept of 
Mayr, 1981), and describe unambiguous 


diagnostic characters (i.e. the phylogenetic 
species concept of Cracraft, 1987). In addition, as 
they are ubiquitous, allozymes offer a yardstick 
to compare levels of evolutionary divergence in 
relation to taxonomic rank in widely different 
taxonomic groups. Through molecular methods, 
it has become easier to verify whether 
ichthyologists, entomologists and spongologists 
infer the same thing when they talk about generic 
taxa in their respective groups. Since 1978, over 
3000 intraspecific and interspecific allozyme 
comparisons have been performed between marine 
populations (literature data based on search on the 
Aquatic Sciences and Fisheries Abstracts 
database, between 1978 and 1998). The most 
commonly used measure of genetic similarity is 
the index of gene identity (7; Nei, 1972), which 
varies from 1.0 (=complete identity) to zero. An 
analysis of the large database of genetic studies, 


mostly for terrestrial vertebrates and Drosophila, 
demonstrated that mean levels of gene identity 
were, as expected, very different when 
conspecific populations, congeneric species or 
confamilial genera were compared (Thorpe, 
1982; Thorpe, 1983). It was shown that less than 
5% of all conspecific comparisons fell below an 
identity level of 0.8 (Thorpe, 1982). 
Consequently, the / valuc of 0.8 has been used as 
a threshold for deciding about specific 
differentiation using allozyme data to define 
species, especially for comparing allopatric popul- 
ations, where the more straightforward usc of 
diagnostic loci (sensu Ayala, 1983) is not possible, 
and the biological species concept (Mayr, 1981) is 
not practical (Aron & Solé-Cava, 1991; Claridge 
etal., 1997). However, that value may be still too 
high for making decisions about the taxonomic 
rank of some marinc invertebrates from 
geographically distant populations. This is 
because the number of allozyme loci detectable in 
marine invertebrates is usually smaller than in 
other organisms, with a consequent increase in 
the variance of estimates of gene identity (Nei, 
1978), and also because gene flow is expected to 
be limited by geographical distance, with a 
consequent lowering of gene identities (Palumbi, 
1992). Considcring that decisions about species’ 
borders in complex groups, using genetic 
attributes, are best taken using what has become 
known as ‘fuzzy logic’ (Van Regenmortel, 
1997), the use of a threshold value becomes very 
important for the comparison of allopatric sponge 
populations. 


Allozyme electrophoresis was first employed 
for molecular systematics of sponge populations 
by Solé-Cava & Thorpe (1986) and recently for 
sponge population genetics (Benzie et al., 1994). 


Molecular data are also very uscful for 
inferring patterns of genetic flow linked to larval 
dispersal (Burton, 1996). Sponge larvae are 
usually short lived (e.g. Borojevic, 1970; Fry, 
1971; Sara & Vacelet, 1973), which suggests that 
geographical distance could determine levels of 
gene differentiation in sponge populations. On 
the other hand, the pattern of gene flow observed 
in many marine invertebrates is often chaotic, 
depending mostly on rare but long-ranging 
broadcasting events (Johnson & Black, 1984). It 
would be interesting, therefore, to verify whether 
gene flow among sponge populations is also 
chaotic or supports the ‘isolation by distance’ 
model of genetic differentiation (Wright, 1978). 


It has been suggested that Porifera might 
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display much higher levels of interspecific gene 
divergence than other invertebrates, possibly due 
to the presence, in the former, of high levels of 
gcne variation (Solé-Cava et al., 1991a; Klautau 
et al., 1994; Boury-Esnault et al.,1999). If this is 
true, then a re-calibration of the threshold value 
of conspecific gene identity should be performed, 
in order to reduce possible type I crrors (i.e. 
deciding that putative species are different when 
they are not), due to a shift in gene identities 
between sponge populations in relation to other 
organisms. This calibration would be 
fundamental both for the analysis of cvolutionary 
rates in the Porifera and for the continuing study 
on putative cosmopolitanism in the group. 


The aims ofthis paperare to: 1) correlate levels 
of intraspecific gcnc identity with geographical 
distance, in order to cstimate the importance of 
larval dispersal to thc composition of sponge 
populations; 2) verify whether patterns of 
interspecific gene similarity in sponges are indeed 
different from those of other marine invertebrates; 
and 3) re-evaluate the threshold genc identity 
value for making taxonomic decisions for sponges. 


MATERIALS AND METHODS 


Data were gathered from the literature and 
from unpublished studies made by our laboratory 
(see references listed in the table legends). 
Whenever necessary, values of mean 
heterozygosity and genetic identity (Nei, 1978) 
were calculated from tables of gene frequency. 
Geographical distanccs were measurcd as the 
shortest distances by sea, using a large scale map 
(1 cm=60km; Christie et al., 1995). The possible 
relationship between pairwise geographical and 
genetic distances for intraspecific populations 
was tested using a Mantel test, with 1,000 
rcplicates (Sokal & Rohlf, 1995). Pooled data of 
pairwise gene identity measures of intraspecific, 
interspecific and intergeneric comparisons were 
used to construct frequency histograms, in a 
similar way as those built by Thorpe (1982, 
1983). The significance of differences between 
mean identity levels in interspecific (intra- 
generic) and intergeneric comparisons was tested 
using a Mann-Whitney U test (Sokal & Rohlf, 
1995). 


RESULTS 


From all available data, 87 intraspecific, 55 
interspecific and 8 intergeneric comparisons were 
compiled (Tables 1-3 respectively). No significant 
correlation (Mantel test; P>0.40) was found 
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TABLE 1. Levels of gene identity between conspecific populations. Key: Km, distance in kilometers; NL, 
number of loci; /, unbiased mean genetic identity (Nei, 1978); H, mean Hardy-Weinberg expected 
heterozygosity (Nei, 1972). References: 1, Benzie et al. (1994); 2, Klautau et al. (in press); 3, Cristiano Lazoski 


(unpublished results); 4, Solé-Cava et al. (1992); 5, Boury-Esnault et al. (1992); 6, 


Boury-Esnault et al. (1999); 8, Sara et al. (1992). 


Bavestrello & Sara (1992); 7, 


Species Locality 1 Locality 2 | Km NL I h Ref 
or MN Willis Island (Aust) Middle Island | &L (| 46 1.00 0.19 1 
Willis Island Magdelaine | 44 | 9 i 0.90 0.26 1 
Middle Island (Aust) Magdelaine 32 8 d 0.89 0.27 1 
Lihou NE (Aust) Lihou SW | 80 ^ 0.93 0.22 d 
| Magdelaine (Aust) Lihou SW 175 A 0.69 0.28 1 
Magdelaine Lihou NE. 9 0.77 0.22 E 
Willis Island Liheu SW 210. a 0.59 025 J 
| Middle Island Lihou SW 210 ^ 0.62 0.21 1 
Willis Island Lihou NE 245 9 0.88 0.20 | 1 
Middle Island Lihou NE 245 9 0.8 os 
2. Chondrilla nucula .. Marseille (Fr) Ligurian (1t) 350 9 0.91 0.11 2 
|3. Chondrilla sp.3 Anjos (Braz) Praia do Forno 2 9 0.69 0.27 2 
| Búzios (Braz) Anjos ja 44 9 0.87 0.56 2 
Buzios Praia do Forno 44 9 0.95 0.56 2 
Itacuruca (Braz) Picinguaba E $ 9 0.95 0.24 2 | 
Búzios Itacuruca — 240 9 0.95 0.30 2 } 
Anjos Itacuruca 240 9 0.92 0.27 2 | 
Praia do Forno (Braz) Itacuruca 240 9 0.93 0.30 2 | 
Picinguaba (Braz) B 1 280 9 0.91 0.22 2 | 
Anjos Picinguaba 300 9 0.95 0.22 2 
Praia do Forno Picinguaba 300 9 | 0.84 025 (| 2 
Büzios | Picinguaba 310 9 0.89 | 0.25 2 
I i Ttacuruca Ilha do Mel 340 9 0.8 0.27 2 | 
Anjos Ilha do Mel 560 9 0.89 0.25 2 
Praia do Forno Ilha do Mel 560 9 0.98 0.27 2 ; 
Búzios Ilha do Mel | 700 Éd 0.93 0.28 2 
IL Noronha (Braz) Búzios 2400 9 0.84 | 027 2 
| Noronha Anjos 2400 9 0.88 0.31 2 
Noronha Praia do Forno 2400 9 2 
Noronha Itacuruca 2600 g 2 
Noronha Picinguaba 2700 9 2 
Noronha Ilha do Mel 3200 9 2 
i pu La Ciota (Fr) Callelongue 17 13 1.00 0.14 3 
1 La Ciota Endoume 25 13 1.00 0.16 3 | 
Callelongue (Fr) Endoume M 8 13 0.89 0.12 _ 3 
La vesse (Fr) Endoume 30 13 0.95 0.11 3 
La vesse La Ciota 3e 13 0.97 0.12 3 
1 = La vesse | Callelongue 2 | Si | un 0:00. u^ "0:08 A 3 
5. Chondrosia sp. Bermudas Recife 6640 13 0.89 0.27 3 
| Bermudas | Büzios 8300 13 0.95 0.33 3 
Bermudas Forno 8330 13 0.95 0.30 3 
Bermudas — |. Angra 8600 13 0.92 0.28 3 
Recife (Braz) . Büzios 1860 13. 0.94 0.27 3 
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TABLE 1. Continued. 


5. Chondrosiu sp. (cont.) Recife — Forno 1890 13 
Recite Angra ul 2160 13 
d um Buzios ic Forno |. 30 13 
Büzios Angra 300 13 
Forno (Braz) Angra 270 13 
6. Collospongia auris Willis Island Middle Island 8.7 6 
Willis Island | Lihou SW | 210 6 
Middle Island Lihou SW L 210 6 
OTI. La Vesse Riou (Fr) 25 16 4 
8. Oscarella lobularis La Vesse Riou 25 16 1.00 0.11 4 
La Vesse Riou dj] 25 12 0.98 0.12 | 5 
9. Petrosia clavata Paraggi (It) Zoagli (It) 2 DUCI EE 0.12 6 
10. Petrosia ficiformis Paraggi Zoagli 2 9 0.90 0.09 6 
W SHOE Willis Island — | — Middle Island 8.7 6 0.96 0.38 1 
Willis Island Holmes 210 6 0.89 0.35 1 
Willis Island Lihou SW 210 i 6 0.85 Il 0.26 l 
i Middle Island Holmes 210 6 0.86 0.34 l 
| Middle Island Lihou SW 210 6 5 1 
Dat Holmes (Aust) Osprey | 300 6 n 1 
ka p Holmes Lihou SW i 370 6 
Willis Island Osprey 430 6 
| Middle Island Osprey 430 6 i 
Osprey (Aust) Lihou SW 630 6 0.49 0.24 | 1 
|. Tiyiiospangia Willis Island Middle Island 8.7 6 0.91 0.30 1 
Diamond (Aust) Lihou SW 50 6 0.89 
] Lihou NE (Aust) Lihou SW 80 6 0.93 
| Diamond Lihou NE 120 | 6 0.86 
IM o Willis Island Diamond 175 6 0.96 
| m ^ Lihou NE Marion 175 6 0.91 
| | — Lihou SW (Aust) Marion 15 | 6 0.99 
j Willis Island Holmes 210 6 0.96 
E Willis Island Lihou SW 210 6 0.85 
= . Middle Island Holmes 210 6 0.90 
E Middle Island Lihou SW 210 6 0.76 
Diamond Marion 220 6 0.93 
z Willis Island Lihou NE 245 6 | 085 
Middle Island Lihou NE 245 6 0.80 
| Holmes Diamond 300 B 6 0.93 
Holmes Lihou SW 370 6 0.85 
| Willis Island Marion 380 a} 6 0.87 
L Middle Island Diamond 380 6 0.91 
Middle Island Marion 380 6 0.80 
Holmes Lihou NE 420 6 0.85 
= ] Holmes Marion 500 6 0.85 
13. Spirastrella hartmani San Blas | (Pan) San Blas 2 l 8 0.95 
San Blas 1 Galeta (Pan) 8 | 0.87 
San Blas 2 Galeta 8 0.95 
14. T. citrina Marsala (1t) Torbay (GB) 3600 11 0.74 
L i Average - 0.89 
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FIG. 1. A, Relationship between geographical distance 
(in logio km) and pairwise gene identities (Mantel test; 
1,000 replicates; P>0.40). B, Frequency histogram of 
gene identity (/) and taxonomic rank for species of 
sponges. C, Frequency histogram of gene identity (/) 
and taxonomic group. Group 1 - Chondrosia, 
Suberites, Petrosia, Plakina and Phyllospongia; 
Group 2 - Axinella, Chondrilla and Clathrina; Group 
3 - Oscarella, Cinachyrella, Tethya, Cliona and 
Spirastrella. 


between geographic distance and genetic identity 
(Fig. 1A). 


The empirical frequency distribution of 
intraspccific (Table 1), interspecific (Table 2) and 
intergeneric (Table 3) gene identities studied on 
the different genera of Demospongiae and 
Calcarea (Fig. 1B) was similar to that found for 
other organisms (Thorpe & Solé-Cava, 1994). 
The average of / over all interspecific sponge 
comparisons was 0.42, which is similar to that 
found among other marine invertebrates 
(720.54). However, the distribution of 
interspecific pairwise gene identities in sponges 
was bimodal (Fig. 1C). Species of some genera 
were consistently highly divergent (/«0.30; 
‘Group 3°: Cinachyrella, Oscarella, Cliona, 
Spirastrella and Tethya), whereas others were 
within the normal range of gene divergence 
(0.40</<0.80; ‘Group 2’: Chondrosia, Suberites, 
Petrosia, Plakina and Phyllospongia). 
Furthermore, in the genera Axinella, Chondrilla 
and Clathrina (‘Group 2’), species displayed both 
low and high levels of genetic differentiation in 
relation to their congeners (0.13</<0.82). Some 
supposedly congeneric species had significantly 
lower (Mann-Whitney U test, z=2.94; P<0.004) 
levels of gene identity (mean /-0.16; Table 2), 
than species of different genera (mean /=0.30; 
Table 3). However, because genetic analyses 
have so far only focused on taxa with depauperate 
morphological characters or other groups 
presenting difficult systematic problems for 
Porifera, a complete pattern cannot be provided 
by the available data. 


DISCUSSION 


Two very interesting results are evident from 
the gene Identity analyses. 1) Generally, levels of 
gene identity were not correlated to geographic 
distance (i.e. it appears that potential for dispersal 
is not a key component in the structuring of 
sponge populations). 2) Levels of interspecific 
gene identities in the few sponge taxa so far 
examined are within the normal range found 
between species of other invertebrates, although 
some sponge genera have species that are 
extremcly divergent from each other. 


The low correlation observed between 
geographical distance and gene identity of intra- 
specific populations suggests that the length of 
larval life is not an essential factor in the 
structuring of sponge populations. Episodic 
recruitment events by rafüng or some forms of 
asexual rcproduction may play a more important 
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TABLE 2. Levels of gene identity between congeneric species. Group 1, Chondrosia, Suberites, Petrosia, 
Plakina and Phyllospongia; Group 2, Axinella, Chondrilla and Clathrina, Group 3, Oscarella, Cinachyrella, 
Cliona, Spirastrella and Tethya. See Table | for key to abbreviations. References: 1, Cristiano Lazoski 
(unpublished results); 2, Bavestrello & Sarà (1992); 3, Benzie et al. (1994); 4, Muricy et al. (1996); 5, Solé-Cava 
& Thorpe (1986); 6, Solé-Cava et al. (1991b); 7, Klautauet al. (in press); 8, Klautau et al. (1994); 9, Solé-Cava et 
al. (19912); 10, Lazoski et al. (in press, this volume); 11, Barbieri et al. (1995); 12, Boury-Esnault et al. (1992); 
13, Solé-Cava et al. (1992); 14, Boury-Esnault et al. (1999); 15, Sarà et al. (1992); 16, Sara et al. (1993). 


|. Group Species | Species 2 NL 1 | Ref 
1 Chondrosia reniformis Chondrosia sp. 12 0.48 I 
1 Petrosia ficiformis Petrosia clavata 9 0.62 3 
i Phyllospongia lamellosa Phyllospongia ulcicornis 6 0.50 3 
1 E Plakina «A Plakina sp.C il 0.49 1 4 
l Plakina A Plakina sp.D ll 0.73 Jie 4 
1 Plakina A Plakina trilopha 11 0.83 4 
1 Plakina C Plakina sp.D il 0.79 4 
1 i __Plakina monolopha Plakina sp.C 11 0.35 4 
1 Plakina monolopha_ Plakina sp. D 1l | 0.58 4 | 
1 Plakina monolopha Plakina trilopha ll 0.61 4 
1 Plakina monolopha Plukina sp.A 11 0.66 4 
1 Plakina trilopha Plakina sp.C ii 0.54 4 
pio Plakina trilopha Plakina sp.D 11 0.61 | 4 | 
1 Suberites pagureorum Suberites luridus 19 0.66 5 
| 1 Suberites pagureorum Suberites rubrus 19 0.67 5 
1 Suberites rubrus Suberites luridus 19 0.98 5 
2 . .. Áxinella damicornis Axinella verrucosa 8 0.13 6 | 
2 Axinella damicornis Axinella sp. 8 0.70 6 
2 Axinella verrucosa Axinella sp. 8 0.13 6 
2 Chondrilla nucula Chondrilla sp.4 (Salvador) 9 0.23 7 
2 Chondrilla nucula Chondrilla sp.1 (Noronha) 9 Mt 0.28 ie 7 
2 Chondrilla nucula l Chondrilla sp.3 (Brazil) 9 0.33 7 
2 Chondrilla nucula Chondrilla sp.2 (Panama) 9 0.53 7 
2 Chondrilla sp.1 (Noronha) Chondrilla sp.2 (Panama) 9 0:32 [ 1 
2 Chondrilla sp.| (Noronha) Chondrilla sp.3 (Brazil) DU O48 j 7 
2 | Chondrilla sp. 1 (Noronha) Chondrilla sp.4 (Salvador) 9 | 0.58 7 
| 2 Chondrilla sp.2 (Panama) Chondrilla sp.3 (Brazil) 9 024 17 
| 2 Chondrilla sp.2 (Panama) Chondrilla sp.4 (Salvador) 9 0.25 7 
| 2 Chondrilla sp.3 (Brazil) Chondrilla sp.4 (Salvador) 9 0.30 7 
| 2 Clathrina aspina Clathrina ascandroides 9 A 0.57 ] 8 
| 2 Clathrina aspina Clathrina cvlindractina 9 0.65 | 8 
|2 Clathrina aspina Clathrina primordialis 9 | 0.82 8 
2 | Clathrina brasiliensis allie Clathrina cylindractina 9 ii 0.43 Al 8 
2 Je Clathrina brasiliensis Clathrina ascandroides _ e 9 0.43 8 
2 Clathrina brasiliensis Clathrina primordialis 9 0.55 8 
2 Clathrina brasiliensis Clathrina aspina 9 0.69 8 
2 Clathrina cerebrum Clathrina brasiliensis 7 0.29 9 
2 Clathrina clathrus Clathrina aurea 11 0.13 9 
AM Clathrina cylindractina Clathrina ascandroides 9 0.43 8 
[2 Clathrina primordialis Clathrina ascandroides 9 0.44 8 
l2 Clathrina priniordialis Clathrina cylindractina 0.65 8 
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TABLE 2. Continued. 


3 Cinachyrella apion Cinachyrella alloclada 19 0.27 si 10 | 
13 Cliona viridis Cliona nigricans 4 | 0.00 _ Elias all. 

3 Oscarella lobularis Oscarella tuberculata 16 Oe 12,13 

3 Spirastrella subogae S. hartmani 8 0.12 a - 

3 Tethya citrina Tethva aurantium 11 0.18 15 

3 Tethya citrina [i Tethya norvegica 11 | 0.20 L 15 

3 Tethya norvegica Tethya aurantium H 0.10 7 AIS 

3 Tethya orphei Tethya robusta "B" 8 0.01 E 16 ] 

3 Tethya robusta "A" Tethya robusta “B” 8 0.18 16 

3 Tethya robusta "A" | — Tethva orphei —— 8 L 0.27 16 

3 Tethva seychellensis Tethya robusta “B”(Red Sea) | 8 | 0.03 I 16 | 

3 Tethya seychellensis Tethya orphei | 8 0.19 16 

3 Tethya sevchellensis Tethya robusta "A" (Maldives) 8 0.28 16 


role in sponges, as observed in other marine 
invertebrates (Johnson & Black, 1984; Johnson 
et al., 1993; Burnett et al., 1995). This indicates 
that sponges follow the islands model, rather than 
the isolation by distance model of genetic 
differentiation (Wright, 1978). Levels of 
population structuring, measured both by 
pairwise gene identities (Nei, 1972) and FST 
inbreeding indices (Wright, 1978) are very high 
in sponge species (Fs7=0.05—0.36; Benzie et al., 
1994; Klautau, unpublished results). This 
indicates that sponge larvae either have low 
capacity for dispersal, or are philopatric, as also 
observed for some species of ascidians (Grosberg 
& Quinn, 1986). In any case, the high levels of 
population structuring observed indicate that 
sponge populations are continuously diverging 
genetically even over small geographic scales. 
Possible consequences of the high level of 
population differentiation are the adaptation of 


TABLE 3. Levels of gene identity between confamilial 
genera. See Table ! for key to abbreviations. 
References: |, Benzie et al. (1994); 2, Solé-Cava et 
al. (1992); 3, Guilherme Muricy & Antonio 
Solé-Cava (unpublished results). 


| Genus 1 Genus 2 NL 1 Ref 

| Phyllospongia | Carterospougia | 6 _ 0.32 g 
Phyllospongia Collospongia 6 0.19 1 
Carterospongia | Collospongia 6 0.20 1 
Oscarella Corticium 16 0.32 2 

j Osearella _ Pseudocorticium | 16 | 0.28 2 
Corticium Pseudocorticium 16 | 0.47 2 | 
Plakina Oscarella 11 0.22 3 
Plakina Corticiun 11 0.30 3 
Plakina Pseudocorticium 11 0.40 3 


local populations to micro-environmental 
conditions, and the scope for a high speciation 
rate in sponges (Benzie et al., 1994). 


In general, the frequency distribution of the 
values of gene identity, in relation to taxonomic 
rank in sponges (Fig. 1 B), shows a similar pattern 
to that observed for other species of animals 
(Thorpe & Solé-Cava, 1994). The main 
differences observed were a slight shift to the left 
in the distribution of intraspecific gene identities, 
and the bimodal distribution of interspecific gene 
identities (Fig. ! B). The higher levels of 
intraspecific differentiation may be related to 
high levels of gene variation (Skibinski & Ward, 
1982) as those usually observed in sponges 
(Solé-Cava & Thorpe, 1989; Solé-Cava & Thorpe, 
1991), although no significant association 
between heterozygosity and gene identity was 
observed for the sponge data (Table 1; Spearman's 
Rank Correlation, P>0.10). The bimodal 
distribution of interspecific gene identities is morc 
puzzling, and seems to result from different 
patterns of gene divergence in different sponge 
genera. The genera analysed can be roughly 
broken into three groups in relation to levels of 
interspecific gene identities: 1) genera whose 
species have similar levels of gene identity as 
other invertebrates (Chondrosia, Petrosia, 
Phyllospongia, Plakina and Suberites); 2) genera 
where some pairwise species comparisons give 
very low identity values (1«0.3), whereas others 
have levels of gene identity comparable to those 
of other organisms (0.4<I<0.8 sensu Thorpe, 
1983; Thorpe & Solé-Cava. 1994) (Axinella, 
Chondrilla and Clathrina); and 3) genera where 
interspecies comparisons consistently give 
extremely low («0.3) identity values 
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[ Corrachyretlla, Cliona. Oscarelfa, Spirastvella and 
Tethva). The ihree groups are significantly 
different from each other (Mann Whitney's U 
test, P<0,001 for each pairwise comparison]. In 
relation to levels of intra and intergeneric 
divergence, group | genera showed a profile 
similar to that observed in other animals, that is, 
species of group 1 were more similar to each other 
than wo species of other generi (Mann- Whilnes 's 
U dest, P«tUL0001). On the other hand, levels of 
gene identity between group 3 species were lower 
ihan those observed between different genera of 
invertebrates (Fig. 1B-C; Thorpe, 1983). More 
interestingly, they were also significantly lower 
(Mann-Whitney's U test, P<0.004) than those 
found between different genera of marine 
sponges (Table 3). 


Consequently, what is the threshold gene 
identity value used for deciding about specific 
differentiation in sponge populations ? 
Taxonomie decisions about the comparison of 
morphs living in sympatry should be based on the 
presence of diagnostic loci (sensu Ayala, 1983; 
ne, loci for which the probability of wrongly 
identifying one individual as belonging to one 
species ofa pair is smaller than 1%), rather than 
using gene identities. The use of diagnostic loci is. 
preferred because it is more consistent from the 
theoretical point of view, and decisions based on 
them amalgamate the power of virtually all 
currently accepted species concepts (Claridge et 
al., 1997). Diagnostic laci can, of course, also be 
found between allopatric populations, but in that 
case making taxonomie decisions about their 
conspecificity is not as straightforward, since 
allopatric populations are expected to diverge if 
levels of gene flow are not very large (Wright, 
1978). Under the phylogenetic species concept 
[Craeraft, 1987), diagnostic loci indicate 
independent evolution, and therefore speciation. 
However, given the very low levels of gene Now 
thal seer La exist between sponge populations, if 
we make an orthodox use of (hat concept we will 
be forced to create new species of sponges for 
almost every allopatric population that we 
analyse. A more reasonable alternative is to use, 
for the comparison of allopatric populations, 
levels of gene identity, since they take into 
account the overall level of divergence, on which 
diagnostic loci do have a heavy weight (/ for 
diagnostic loci = 0). but that is less biased hy 
Spire events of selection or drift, Given the 
shift to the left in the intrageneric gene identity 
distribution, a value of about 0.7 of gene identity 
can be chosen as a threshold for making decisions 


MEMOIRS OF THE QUEENSLAND MUSEUM 


about interspecitic dilferentiation berween 
allopatric sponge populations, We chose this 
value because it corresponds to the point where 
the distribution curves of intraspecific and 
interspecific gene identities meet, clearly 
separating the two groups (Fig. ! B), If we observe 
the distribution of interspecific identity levels 
(Tahle 21 and consider only the comparison of 
sympatric populations, thus avoiding the 
potential circularity of using allopatric 
comparisons, we can see that the average of 7 is 
0,52, and only 16% of the interspecilic gene 
identity values are above 0.70. Using this 
threshold value, Phyllospengia aleicornis trom 
Lihou and Osprey Reefs (Coral Sea, Australia), 
and some of the populations of Carteriospongia 


fahellifera from Lihou Reef, considered hy 


Benzie el af, (1994) to be conspecific with those 
of Middle Island and Willis Island, would have ra 
he considered as separate species (/=0.49-0,67), 

The belief that sponge species have a higher 
level of genetic differentiation than other 
organisms (Sole-Cava et al, 199 La; Klautau et 
al., 1994; Boury-Esnault etal., 1999) may simply 
he the result of an over-representation of species 
of group 3 in the literature, The picture that 
emerges [rom this study using a larger set oldala. 
indicates that levels of gene divergence among, 
presently recognised sponge genera vary 
broadly, which may be the result of two different, 
but not mutually exclusive, phenomena. || Those 
sponge genera with genetic identities helow 0.3 
are so old that there has been a saturation of eene 
divergence leading to the accumulation of 
homuplasic changes (as discussed by Thorpe. 
1989). 2) Some sponge genera, notably those of 
groups 2 and 3 above, are polyphyletic. 1n the 
first case, allozymes would be considered to he of 
little use above the species level in sponges, but it 
would remain fo be explained why some sponge 
genera can diverge at so different rates 
(Solé-Cava & Thorpe, 1994), In the second 
instance, Some sponge genera require revision, 
and possibly splitting up into smaller, 
monophyletic units. In the first case the 
vongenerie species of group 3 should be much 
morc different from cach other than those of 
different genera of practically all other groups of 
animals (including sponges). 

The high levels of gene divergence observed 
between conspecific sponge populations und 
between species in some genera of sponges should 
be further investigated, as they have important 
consequences for the taxonomic framework for 
the whole group. If the gene identity found 
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between species of Cliona is zero (Barbieri et al., 
1995), and between species of Spirastrella is 
0.12, what would be the identity between Cliona 
and Spirastrella species? Likewise, what levels 
of gene identity would be observed between 
species of Tethya and Tectitethya or Timea ? At 
this very low level of gene identity, intergeneric 
species may have some alleles in common by 
simple homoplasic convergence, due to the 
saturation of possible alleles detectable by the 
technique (Thorpe & Sole-Cava, 1994). Those 
convergent allelcs are often found in taxonomic 
comparisons above the genus level, but their 
presence is usually detected because they conflict 
with a much larger number of true synapo- 
morphies within each genus (Hillis et al., 1996). 
However, given the very low gene identity found 
between spccies of group 3 genera, these few 
convergent alleles could be misinterpreted as 
synapomorphies, and lcad to wrong taxonomic 
conclusions. For example, considering the lack 
of synapomorphies in the molecular data within 
Cliona or Cinachyrella, and the possible alleles 
in common between species from those genera, 
what should be our decision about their 
taxonomic status? Further genetic studies, 
possibly linked to independent DNA analyses, 
are nccded to determine whether allozyme data 
are sufficently objective to distinguish sponges at 
the genus level. 
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